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Abstract: Motivated by the recent proliferation of observed astrophysical anomalies, 
Arkani-Hamed et al. have proposed a model in which dark matter is charged under a non- 
abelian "dark" gauge symmetry that is broken at ~ 1 GeV. In this paper, we present a survey 
of concrete models realizing such a scenario, followed by a largely model-independent study 
of collider phenomenology relevant to the Tevatron and the LHC. We address some model 
building issues that are easily surmounted to accommodate the astrophysics. While SUSY is 
not necessary, we argue that it is theoretically well-motivated because the GeV scale is auto- 
matically generated. Specifically, we propose a novel mechanism by which mixed D-terms in 
the dark sector induce either SUSY breaking or a super-Higgs mechanism precisely at a GeV. 
Furthermore, we elaborate on the original proposal of Arkani-Hamed et al. in which the dark 
matter acts as a messenger of gauge mediation to the dark sector. In our collider analysis we 
present cross-sections for dominant production channels and lifetime estimates for primary 
decay modes. We find that dark gauge bosons can be produced at the Tevatron and the 
LHC, either through a process analogous to prompt photon production or through a rare Z 
decay channel. Dark gauge bosons will decay back to the SM via "lepton jets" which typically 
contain > 2 and as many as 8 leptons, significantly improving their discovery potential. Since 
SUSY decays from the MSSM will eventually cascade down to these lepton jets, the discovery 
potential for direct electroweak-ino production may also be improved. Exploiting the unique 
kinematics, we find that it is possible to reconstruct the mass of the MSSM LSP. We also 
present several non-SUSY and SUSY decay channels that have displaced vertices and lead to 
multiple leptons with partially correlated impact parameters. 



1. Introduction 



Several intriguing observational results from high-energy astrophysics have motivated an ex- 
citing new proposal [1] in which a WIMP-like dark matter (DM) particle at 500-800 GeV 
annihilates primarily into leptons and is charged under a new "dark" force carrier. ATIC [2] 
detects an abundance of cosmic ray electrons between 300 — 800 GeV, while PAMELA [3] 
sees an excess of positrons (but not anti-protons [4]) at 10-100 GeV. Together with the CMB 
haze [5, 6, 7], these observations paint a consistent picture whereby DM annihilates primarily 
into muons and/or electrons [8]. 

There are two sources of tension between these results and more conventional models 
of WIMP dark matter. First, assuming thermal freeze-out, the standard relic abundance 
calculation implies an annihilation cross-section that is at least a hundred times too small to 
explain the lepton excesses observed in astrophysical experiments. A "boost factor," typically 
attributed to local over-densities of dark matter, is often evoked in this case. A second 
difficulty is the non-observation of corresponding excesses in anti-protons [4] and gamma rays 
[9], which puts strong bounds on hadronic channels that are present in many dark matter 
models. 

Motivated by the above considerations, the authors of Ref. [1] outline a scenario in 
which these apparent contradictions are reconciled. They introduce a 500-800 GeV WIMP 
that couples to a GeV scale dark gauge boson that kinetically mixes with the photon of 
the Standard Model (SM) [10] 1 (see Ref. [12] for another recent suggestion with similar 
ingredients). A schematic illustration of this scenario is presented in Fig. 1. The ATIC 
and PAMELA data are explained by DM annihilation into the dark gauge boson which 
subsequently decays into electrons and muons. Elegantly enough, the 0(1) GeV scale plays 
two independent roles. First, the new dark force carrier at < GeV introduces a Sommerfeld 
enhancement [13, 14, 15, 16], giving a boost factor of the right size to enhance the DM 
annihilation cross-section 2 . Second, the absence of anti-protons in the PAMELA observations 
is now simply a result of kinematics [18]. 

The gauge group, Gdark; is a priori unspecified. However, it was observed in Ref. [1] that a 
non-abelian Gdark nicely accommodates the excited dark matter (XDM) [19] and inelastic dark 
matter (iDM) [20] mechanisms. XDM was proposed in order to explain the INTEGRAL [21] 
measurement of the 511 keV gamma-ray line at the center of the galaxy. The iDM scenario can 
accommodate the DAMA /LIBRA measurement of WIMP-nuclei scattering with other direct 
detection experiments [22, 23, 24]. Both XDM and iDM are non-standard WIMP scenarios 
in which a DM ground state can transition to and from new excited states via the emission 
of some field that couples back to the SM. If the DM lives in a multiplet of a non-abelian 
Gdark; then these ground and excited states can be the components of this multiplet, and 
transitions will emit dark gauge bosons that couple weakly to the SM. Independent of the 

1 Ref. [11] analyzes particle physics bounds on such a light vector field and its possible connection to the 
HyperCP anomaly. 

2 See Ref. [17] for an alternative, but related way for producing a large boost factor. 
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Figure 1: A schematic illustration of the minimal setup we consider in this paper. The dark sector 
and the SM are connected through kinetic mixing term suppressed by e < 10 -3 . The dark matter 
multiplet may or may not couple directly to the SM. Supersymmetric extensions of this scenario are 
also discussed. 

results from INTEGRAL and DAMA, we find the possibility of a non-abelian dark sector to 
be intriguing in its own right, with direct implications for the collider phenomenology. Thus 
throughout this paper we consider a dark sector with a non-abelian gauge symmetry that is 
completely broken by some dark Higgs sector 3 . 

In Section 2, we construct a catalog of explicit minimal models. Since Gdark needs to 
include a U(l) factor for kinetic mixing with SM hypercharge, we take Gdark = SU(2) x £7(1). 
Our models differ only in their dark Higgs sectors, which are constructed to break Gdark 
completely and induce all the necessary couplings between the different states of the DM 
multiplet. 

In Section 3, we discuss the mass splittings between the dark matter states. In order to 
obtain the small mass splittings needed for XDM and iDM, we consider DM that is a doublet 
or a triplet under 5C/(2)d ar k- The splittings may be generated radiatively from dark gauge 
boson loops. Another possibility is to generate them through higher-dimensional couplings 
between the dark matter and a single dark Higgs. 

In Section 4, we consider the addition of SUSY to the dark sector. We observe that the 
minimal assumption of kinetic mixing between dark sector and SM hypercharge generates 
an effective FI term in the dark sector that is naturally of the desired scale, O(GeV). This 
term can break SUSY, or even more interestingly can generate a super-Higgs mechanism that 
leaves a supersymmetric dark sector with a ~ 1 GeV gap. Both of these scenarios typically 
result in light fermions that may have an influence on collider physics. We emphasize that 
this is a leading contribution which must be included in any SUSY scenario that includes 
kinetic mixing. Furthermore, within this scheme the DM can easily be a SM singlet, and 
so DM annihilations do not produce SM bosons that would dangerously decay to anti- 
protons that have not been observed by PAMELA. We also investigate the gauge mediation 
scenario originally proposed in Ref. [28] where DM is charged under the SM gauge group. 
An additional complication we address arises because SUSY restricts the form of the scalar 

3 There are strong astrophysical constraints on a long range interaction from unbroken gauge symmetry 
with an unsuppressed coupling [25, 26, 27]. 
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potential which is responsible for breaking Gdark completely. We provide several examples to 
overcome this difficulty. 

In Section 5, we present several benchmark models for the dark Higgs sector. The re- 
sulting spectra of light vector bosons and scalars are explicitly computed and the relevant 
couplings are discussed. 

In Section 6, we investigate the collider signatures of these models 4 . The kinetic mixing 
is the essential gateway to produce and observe dark sector states. Dark gauge bosons can 
be produced in processes analogous to prompt photon production in the SM. They can also 
be produced through rare Z decays. The dark sector states themselves dominantly decay 
into multiple and ft* 2 , which are highly collimated and dubbed "lepton jets" [28]. Due to 
the non-abelian structure of Gdark; these lepton jets typically contain more than two leptons 
each. We discuss the observability of such signals at the LHC and at the Tevatron. We 
find that the cascade decays in the dark sector may result in displaced vertices or possible 
correlations between the fix and the lepton jet. Several such displaced vertices will produce 
uncorrelated impact parameters of decay products. In the case where the dark sector is 
supersymmetric, then it may be possible to detect direct electroweak gaugino production 
with enhanced reach both at the Tevatron and at the LHC. As a bonus, we find that we can 
exploit these cascades to perform absolute mass measurements of MSSM gauginos. Section 7 
contains our conclusions. 

Finally, let us briefly comment on the notational conventions used in this paper. In 
general, symbols referring to elements of the SM will be capitalized — so for example the SM 
hypercharge gauge coupling, gauge field, and field strength will be denoted by gy, B„ and 
B^ u . In contrast, lowercase symbols will refer to elements of the dark sector, so the dark 
sector gauge coupling, gauge field, and field strength will be denoted by g y , and b^ v . We 
will use h, or h! to denote dark Higgses. We denote the dark matter states by ^ and we 
denote the SM and dark photon by 7 and 7', respectively. 

2. The Dark Sector and Symmetry Breaking 

Let us begin by discussing the basic structure of the dark sector models that we will consider 
in this paper. We take the DM to be the lightest (and stable) component of some multiplet 
of the non-abelian group Gdark- As we will discuss in Section 3, such a multiplet is necessary 
if we wish to explain the INTEGRAL and/or DAMA signals along the lines of the XDM and 
iDM proposals of [19, 20, 29]. 

Furthermore, we follow the proposal of [1] in which the SM is coupled to the dark sector 
via a kinetic mixing term between SM hypercharge and a dark sector U (1) gauge field: 

£gauge = -\B^B^ - \b^ V + ^B^ (2.1) 

where B^ v and b^ v are the SM and dark sector hypercharge field strengths, respectively. 
Because this marginal operator preserves all of the symmetries of the SM, it is relatively 
4 We leave any precise matching to astrophysical observations for future work. 
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unconstrained phenomenologically. For a detailed analysis of kinetic mixing and the couplings 
it induces between SM and dark sector fields, see appendix A. 

Since Gdark must contain a U(l) factor 5 , the minimal choice is of COUrse Gdark — SU (2) x 
U(l). Furthermore, if Gdark is broken completely at a scale of ~ GeV, then the resulting 
mass gap will relieve constraints from BBN on the number of relativistic degrees of freedom. 
However, in order to fully break charge, it is necessary to appropriately engineer a dark Higgs 
sector. As we shall see shortly, these scalars must also break a custodial SU(2) in order to 
be phenomenologically viable. The necessity of breaking these symmetries demands a fairly 
elaborate dark Higgs sector. 

First, let us consider the issue of charge breaking. Even for the simplest two Higgs 
doublet model, the criterion for charge breaking is quite complicated [30], for theories with 
more exotic Higgs representations, the space of charge breaking vacua is not even known. In 
appendix B, we present a straightforward method for deriving necessary conditions for charge 
breaking in two higgs doublet sectors, which we applied in order to obtain viable dark sector 
benchmark models. 

Now let us explain the problem of the custodial symmetry. In the spirit of [1], we will 
assume that the DM is a multiplet of Gdark whose components are split in mass. The resulting 
excited and ground states have transitions mediated by dark gauge bosons that need to couple 
to the SM electric current if they are to realize the XDM and/or iDM scenarios (see Section 3). 
However, this mixing can be forbidden by a custodial symmetry of the dark Higgs sector. To 
see why this is the case, consider a model of two scalar doublets. We define w a ^ and b^ to be 
the gauge bosons of Gdark> where is the abelian field which mixes with the SM hypercharge, 
B^. Assuming arbitrary vevs for the scalars, Gdark is broken, and in the {w\, W2, W3, b} basis, 
the gauge boson mass matrix takes the form 



As a consequence of the custodial symmetry present in any theory of only scalar doublets, the 
diagonal entries Wi are all equal. Applying a custodial SU (2) transformation, we can rotate 
the components Aj completely into the direction. This yields a mass matrix which has 
a manifest £7(1) symmetry that acts as a phase rotation on w± = w\ ± iu)2 (note that the 
gauged "electromagnetism" can still be broken while preserving this U(l)). Under this U(l) 
the components of the DM multiplet have distinct charges — consequently the gauge bosons 
that mediate transitions among these states must also be charged, so they can only be the 

5 It is actually possible to achieve mixing without an abelian factor through an S parameter type operator 
Tr B M ", where $ is some operator that transforms as an adjoint of the non-abelian group. In this 

paper we keep the abelian factor in order to investigate the collider signatures of the more general gauge 
group structure and ignore the existence of such operators. That is certainly justified in the case where no 
fundamental adjoints are present and the contribution is subleading. 
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w±. However, w± have no components in the b direction, so they do not kinetically mix with 
SM hypercharge and thus cannot decay to SM particles. 

Because the custodial symmetry is broken explicitly by the dark hypercharge, the cou- 
plings to the SM that are excluded at tree-level by this symmetry will be generated at one 
loop. Indeed, this may actually be desirable, since it generates an effective coupling for the 
DAMA transition that is suppressed beyond the e 2 from the kinetic mixing. Another possi- 
bility, considered below, is to include additional Higgses that break custodial symmetry at 
tree-level. 

In the case of SUSY models, we will be forced to significantly enlarge the Higgs sector. 
This is because many of the difficulties that arise in the non-SUSY case are exacerbated with 
the additional constraints imposed by SUSY. Moreover, in SUSY, all scalars are complex, 
which forces us to promote real Higgs triplets to complex Higgs triplets. This, along with 
the constraint of anomaly cancellation implies somewhat of a proliferation of Higgses in these 
theories. 

In what follows, we enumerate several types of scalar sectors that break dark charge as 
well as custodial symmetry. We focus on models with the intention of later extending them 
with supersymmetry. 

2.1 Doublet Models 

A theory of one Higgs doublet is incapable of breaking charge, so we consider two doublets 
hi and h 2 with quantum numbers 2_ 1 / 2 and 2 l j 2 under Gdark = SU(2) x U{\). A general 
renormalizable scalar potential that breaks charge is given by 6 , 

V(h u h 2 ) = y (M 2 - M 2 ) 2 + y (N 2 - H 2 ) 2 (2.3) 
+ A 4 \hjeh 2 - viv 2 cosa\ 2 + A3 (\hi\ 2 - \vi\ 2 ) (\h 2 \ 2 - | ^2 1 2 ) 

with, 

Ai > 0, A 2 > 0, ^f>^\2 + A 3 > 0, A 4 > (2.4) 

v 2 is complex and charge is broken when < | cos a\ < 1. 

In the MSSM the conditions of Eq. (2.4) are violated at tree level. From the D-term 
contributions to the scalar potential we have X\ = X 2 = —A3 = (g 2 + g' 2 )/8. The inequality is 
saturated and the potential in Eq. (2.3) degenerates and contains a flat direction. To avoid 
such flat directions in the MSSM one must usually evoke a condition on the quadratic terms 
in the potential. Such potentials cannot be placed in the form of Eq. (2.3) and charge is not 
broken. Therefore the usual super symmetric two doublets model will not suffice and we need 
additional contributions to the scalar potential in order to satisfy the condition, Eq. (2.4). 

6 This is not the most general renormalizable scalar potential possible. One can add three more terms, 
|fti| 2 |ft2| 2 , \hi | 2 /ife/i2, and \h2\ 2 hfeh2 which are consistent with all the symmetries. However, these simply 
complicate the potential and are not required for breaking charge. A more general analysis of the vacuum 
structure can be found in Ref. [30] 
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In addition, since this model has a custodial symmetry, it fails to have proper mixings 
between the gauge bosons. Nonetheless, since the custodial symmetry is broken by dark 
hypercharge, the gauge boson mixing receives one-loop radiative corrections that break the 
custodial symmetry. From this point of view there is also no reason not to include higher 
dimension custodial violating operators that can be generated if heavy (triplet) states have 
been integrated out. In fact, we include such irrelevant operators in the benchmark model of 
Section 5.1.1. 

Here we also note the presence of an unfortunate Z2 symmetry that is present in the 
tan/3 = 1 limit. This symmetry needs to be broken since it prevents two of the dark gauge 
bosons from coupling to SM electric charge (see appendix C). 

2.2 Doublet/Triplet Models 

An obvious way to break custodial ST/ (2) at tree-level is to augment the two doublet model 
with a light triplet of SU{2). For instance, consider a model of one doublet, h, and one 
real triplet, <I>, with dark quantum numbers 2 1 / 2 and 3o, respectively. In order to realize a 
charge breaking angle between the doublet and triplet, we include the following two operators: 
h^&h and h T e&h. Since the latter has nonzero hypercharge, we must multiply it by a new 
hypercharged singlet, S, in order to include it in the potential: 

V(h, S) = ^ (|/f - M 2 ) 2 + ^ (Tr [<M>] - M 2 ) 2 + ^ (|S| 2 - M 2 ) 2 (2.5) 
+ Cl h^h + {c 2 Sh T e<S>h + h.c.) 

Alternatively, we might consider a model with two doublets and one triplet. This is more 
natural if we wish to eventually include SUSY. The scalar potential takes the form: 

V(h u h 2 ,$) = V(h u h 2 ) + ^ (Tr [$$] - |^| 2 ) 2 (2.6) 
+ ah{®hi + c 2 h\$h 2 + (c 3 hj e®h 2 + h.c.) 

where V(h\, h 2 ) is the contribution from doublets alone defined in eq 2.3. 

We can impose an additional % 2 symmetry <I> — ► — <£, that forbids tree- level couplings 
between the triplet and doublets: c\ = c 2 = C3 = 0. This enhanced global symmetry implies 
the existence of two pseudo-Goldstone bosons which obtain masses at one-loop ~ 10 MeV. 
These pseudo-Goldstone bosons will be produced at the bottom of dark sector cascades. 
They decay into leptons through either two off-shell dark gauge bosons or at one-loop (see 
Fig. 10 and the discussion in Section 6.1.1). Either way, the long lifetime causes the pseudo- 
Goldstone boson to escape the detector at colliders. Since those are pseudo-scalars they will 
not contribute to the Sommerfeld enhancement of DM annihilations in the early universe and 
their mass is therefore not bounded by the limits derived in Ref. [31]. 
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3. Dark Matter Mass Splitting 



The authors of [1] observed that a DM multiplet of some non-abelian Gdark; given appropriate 
mass splittings, can in principle realize the XDM explanation of INTEGRAL [19] and also the 
iDM mechanism for reconciling the DAMA annual modulation with the null result of other 
direct detection experiments [20, 29]. In this section we briefly review these proposals, and 
discuss concrete ways of generating the appropriate mass splittings within concrete theories. 

The INTEGRAL collaboration has provided an extremely refined measurement of the 
511 keV line of positronium annihilation coming from the galactic center. In the XDM 
scenario, WIMPs in the galactic center scatter into an excited state, lying ~ 1 MeV above 
the ground state. The excited state then de-excites into e + e~ which provides the excess 
positrons needed. In terms of model-building we need a splitting of ~ 1 MeV between two 
states in the DM multiplet. Transitions between these two states are mediated by a dark 
gauge boson with some component of the dark hypercharge (which in turn couples to SM 
leptons). 

In contrast, DAMA is a direct detection experiment which seeks to measure the scattering 
of galactic WIMPS off of Nal(Tl). Assuming a standard WIMP with elastic scattering, several 
other experiment such as CDMS[22, 23], XENON [24], and ZEPLIN [32] exclude DAMA's 
measured annual modulation by many orders of magnitude. The iDM proposal reconciles 
these experiments by proposing that the WIMP can only scatter off of nuclei through an 
inelastic process by which the the DM is converted into a slightly excited state. Since the 
WIMP kinetic energy is fixed and the threshold for the inelastic transition is dependent on 
the atomic number of the nuclei, this iDM scenario can simultaneously predict a null result 
at CDMS and a positive result at DAMA 7 . Considering fermionic DM, this scenario can be 
accommodated by including a mass splitting of around ~ 100-150 keV [20, 29] between the 
lightest two Majorana states of the fermion. The bottom line for model building is that 
to evade CDMS and CRESST[33] bounds, the DM must be split from the next heaviest 
Majorana state by at at least 100 keV. 

Before we consider mechanisms for generating the required splittings, we must ascertain 
that there is no elastic scattering which would have been seen in direct detection experiments. 
One possibility is to begin with Majorana dark matter in a real representation of the dark 
gauge symmetry. Gauge bosons then couple different states of the multiplet and radiative 
corrections, to be discussed in Section 3.1, can split the masses of these states. But if the dark 
matter begins in a complex representation, for example if it has dark or SM U(l) charge, then 
it must be Dirac-like at high-energies. Then the model is already excluded by direct detection 
experiments since the elastic scattering of Dirac-like dark matter is not sufficiently suppressed 
unless e < 10~ 6 . However, it is possible to split the masses of the Majorana components of 
the Dirac fermions by using the same scalar sector that is responsible for breaking dark gauge 
symmetry. For instance, if we imagine that cj> is some scalar singlet whose vev breaks global 

7 XENON and ZEPLIN, which both use Xe as a target should be able to exclude the iDM scenario, but at 
the moment these experiments are background limited [29] 
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fermion number and U(l) y , we can add a term such as, 

Avlajora„a = <^ + 0** C * C (3.1) 

where ^ and ^ c are the Weyl components of some DM multiplet. If <j> develops a vev of order 
~ GeV, it will generate a Majorana mass splitting that forbids any elastic scattering and 
evades direct detection bounds. 

Another possibility is to use a higher dimensional operator with a dark sector doublet, 

h, 

^Majorana = —h^h, (3.2) 

Mx 

where Mx ~ TeV. In this case, the Majorana splitting is of order ~ MeV which again 
kinematically forbids elastic scattering. If DM is charged under the SM as part of a 5 + 5 
multiplet, then a dimension 6 operator is required to contract both dark hypercharge and 
SM quantum numbers. For example, we can use the operator 1/M X H^fyHcf), where H is 
the SM Higgs, is a singlet that soaks up \I/'s dark hypercharge and gets a vev at ~ 1 GeV, 
and Mx ~ TeV. This leads to Majorana splitting of order ~ GeV, which forbids elastic 
scattering. 

Any of the possibilities mentioned above can be employed to evade direct detection from 
CDMS. In the next subsection, we consider two possible way for generating the appropriate 
~ 1 MeV and ~ 100 keV mass splittings necessary for XDM and iDM. 



3.1 Radiative Splitting 

As is well-known [34], spontaneous symmetry breaking of a non-abelian gauge group generates 
radiative mass splittings within a multiplet of the symmetry. We take the DM multiplet to 
have mass ~ 500 — 800 GeV, and to be charged under the dark SU(2) x J7(l). As discussed in 
section 2, realistic dark sectors must break charge and custodial SU(2), but to develop some 
intuition about the radiative mass splittings, we will begin by considering the limit where 
these symmetries are preserved. In this limit the mass splittings among the multiplet take a 
particularly simple form, 

„,dark 

Arriij = -—( q *-qj)M g (3.3) 

dark 

- \~ m 2 - (Tj) 2 ) (M z - M w ) , 

where we define a2 ark , and a dark as usual with respect to SU{2) x U(l) couplings. The charges 
are qi = Tf + Y and Tf is the i th eigenvalue of the third ST/ (2) generator. In the more general 
limit where charge and custodial symmetry are broken, one must use the appropriate vector 
boson mass eigenstates and their couplings to the fermions in order to compute the mass 
correction (Eqs. D-l). This is a straightforward computation, however, in general it does not 
yield a simple analytic result. Nevertheless, it is clear that there are two factors which control 
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Figure 2: The ratio of the XDM splitting to the iDM splitting as a function of triplet dark matter 
U(l) y hypercharge. The green horizontal line indicates the minimum ratio for simultaneously achieving 
both splittings. Red (line) is an example of two Higgs doublets with charge preserved, blue (dashed) 
represents two Higgs doublets with charge broken, and black (dots) adds a Higgs triplet to the previous 
case. For this example, the gauge couplings are g — 0.97 and g y — 0.26, and in terms of Eq. 2.3 we 
have for all three models V\ = 0.9 GeV, V2 = 1.1 GeV and Ai 2,3,4 = 1- Red (line) and black (dots) add 
charge breaking with cos a = 0.75, and for black (dots), in terms of Eq. 2.6, A$ ~ 1, t>$ = 1 GeV and 
the triplet is decoupled from the doublets at tree-level by imposing the discrete symmetry: <f> — > — «&. 

the mass splitting: first, the differences between masses of the vector bosons; and second, the 
couplings of the different members of the representation to the vector bosons. 

As a simple example with all the required splittings and couplings we can consider a 
triplet with hypercharge y = 1/2 — S. We generate both large, AM ~ a da,rk M z and small 
Am ~ £)ev dark M 2 splittings. The correct couplings to account for the XDM and iDM scenarios 
are induced when charge and custodial breaking corrections are included. A realistic model 
certainly need not be based on such odd charge assignments, however, this example serves to 
illustrate how straightforward it is to obtain the correct splittings and couplings. In Figs. 2 
and 3, we consider some of the more general models of Section 2, which include charge and 
custodial symmetry breaking, and we plot the exact ratio of the two splittings relevant to 
XDM and iDM as a function of the parameters. The corrections induced in supersymmetric 
models are discussed in appendix D. 

3.2 Mass splitting from higher dimensional operator 

It is also possible to generate the INTEGRAL and DAMA mass splittings from higher di- 
mension operators alone. The key observation is that dm ~ Aj ark /Mx ~ MeV, which is of 
the desired range. 

As an example, we consider two Weyl fermions Vl/,\E rC which are 2 1 / 2 and 2_ 1 / 2 under 
^dark- It is possible to achieve all the required splittings and transitions with a single scalar 
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Figure 3: Two contour plots of the ratio of the XDM splitting to the iDM splitting for triplet dark 
matter with two Higgs doublets and one Higgs triplet. The shaded regions represent splitting ratios 
where XDM and iDM can be achieved simultaneously. In both plots, the horizontal axis is the dark 
matter U (l) y hypercharge. The vertical axis of the left plot represents the ratio of the triplet to doublet 
VEVs, v$/v, where v 2 = v 2 , + v\ and ($) = w$T 3 . The vertical axis of the right plot represents the 
ratio of dark hypercharge and SU{2) couplings, g'/g. For both plots, the triplet is decoupled from 
the doublets at tree-level by imposing the discrete symmetry: $ — > — $, and in terms of Eq. 2.3 we 
have v\ = 0.9 GeV, V2 — 1.1 GeV, cos a = 0.9 and Ai.,2,3,4 = 1- For the left plot, the gauge couplings 
are g = 0.97 and g y = 0.26. For the right plot, we have also chosen, in terms of Eq. 2.6, A$ = 1 and 
v* = 1 GeV. 



doublet, 

c d M 9 m c + —wi^h + — n c n c + —m c h c ^h + h. c . (3.4) 

M x M x M x 

with \& = (ip u ,ip e ) and *$>h = ipiCijhj and h\ = eijh*-. Once the scalar doublet gets a vev, 
(h) = (0, v), the "neutrino" components of ^> and \I/ C mix through the following matrix, 

V My + X 3 v X 2 v J 

where v = v 2 /Mx ~ MeV. In the limit where Ai = A2 = the states are maximally mixed, 
t/j = ± ipl) /v2, and form a Dirac pair of mass My + \%v which provides the XDM 
splitting when compared with the ip e ,ipl states of mass My. 
With non-zero Ai and A2 we have, 

V'i = cos + sin Otfj- mi = My + |(2A 3 + Ai + A 2 ) (3.6) 

ip' 2 = -sin6>^ + + cos 6ip~ m 2 = -My - |(2A 3 - Ai - A 2 ) 

with sin (9 ~ (Ai — \ 2 )v/4My. The mass difference between the two states is |Ami2| = 
(Ai + \ 2 )v. So, by tuning Ai against X 2 we can achieve Ami 2 ~ 0.1 MeV = 100 keV as 
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Figure 4: The resulting spectrum for a Dirac doublet with majoron coupling. 

required by iDM. The coupling of the mass eigenstates to the dark gauge boson is given by, 

g v VlV - g y V e l* c = g y sine cos eWJ^-fiJrPi) 

- g y cos 2 Oip'Jip' 2 + h.c. (3.7) 

In this case the ratio of the elastic to inelastic coupling is approximately sin# = (Ai — 
\2)(v/M#) ~ 10~ 7 , which is sufficiently suppressed. The spectrum relevant for this case is 
shown in Fig. 4. 



4. Generation of the Dark Sector Mass Scale 

As noted in [28] , a particularly nice feature of a SUSY dark sector is that the GeV scale is 
naturally generated by gauge mediated SUSY breaking from the SM. In this section, we elab- 
orate on this scenario in detail. Furthermore, we propose an even more minimal alternative 
in which "kinetic mixing mediation" breaks SUSY or induces a super-Higgs mechanism at a 
scale of several GeV in the dark sector. As we will discuss, these theories typically have light 
fermions which affect the collider physics. 

For gauge mediation, dark matter itself can act as the messenger if we take it to be 
charged under the SM as part of a 5 + 5 multiplet. Dark matter annihilations then also 
produce SM electroweak gauge bosons, resulting in hadronic channels. But since the GeV 
scale can be generated by kinetic mixing mediation alone, there is no need to charge dark 
matter under the SM. 

Although we focus on kinetic mixing mediation and gauge mediation for the rest of this 
section, and when we construct benchmarks in Section 5, there are other ways to break 
SUSY in the dark sector. We would like to stress that the rest of our paper, in particular 
the model-independent discussion of collider signatures in Section 6, does not depend on 
how SUSY is broken in the dark sector. One alternative is that there is high-scale gauge 
mediation and a GeV scale gravitino [28]. Then SUSY is broken in the dark sector at the 
GeV scale by a "Planck slop." Another possibility is that the dark matter mass is related to 
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the mechanism that sets the MSSM fi parameter, for example due to a superpotential of the 
form: A SH u Hd + A' S^f^> c . A vev for F$ is communicated to the dark sector through gauge 
mediation with dark matter as messengers. 

4.1 SUSY Breaking from Kinetic Mixing 

In [35] it was observed that mixing of gauge boson kinetic terms will induce mixed D-term 
contributions to the action that can communicate SUSY breaking between two sectors that 
are otherwise decoupled. 8 There the authors noted new, possibly dangerous contributions 
to SUSY breaking to the MSSM from this effect. In this section, we use this effect to our 
advantage in order to mediate SUSY breaking from the SM to the dark sector. We should 
note that while we can choose to make this the dominant mechanism for breaking SUSY in 
the dark sector, it is always present at the GeV scale. 9 

As was originally proposed in [1], we have been assuming that the dark sector and the 
SM are coupled via a marginal gauge kinetic mixing between the dark hypercharge, U(l) y 
and the SM hypercharge, U(l)y. If both ?7(l)'s are fundamental, then the kinetic mixing is 
a UV boundary condition, sensitive to physics at the highest scales. Instead, if either U(l) is 
embedded in a GUT, then the kinetic mixing is only induced below the scale of GUT breaking 
by integrating out fields charged under both C/(l)'s. In this case we can estimate its size. In 
particular, heavy fields charged under both the SM and the dark sector will induce a gauge 
kinetic mixing: 

£gau g c = \ j d 2 ° (WyWy + W y W y - 2eW Y W y + h.c.) 
16vr 2 8 \M> 2 / 

where g y and gy are the gauge couplings for the dark and SM hypercharges, respectively, 
M and M' are the masses of components of the heavy particle multiplet. Assuming that 
these mass scales are not too separated and that the gauge couplings are of reasonable size, 
this gives an estimate of e ~ 1CU 3 — 10~ 4 . Interestingly, this not only gives the right scale 
to explain the DAMA cross-section, but also generates a scale of around a GeV in the dark 
sector. Eq. (4.1), along with the Kahler potential, implies a D-term potential: 

U gaugc = X -D\ + X -Dl - eD Y D y + g Y D Y £ Q t \H t \ 2 + g y D y ^ Qi \h\ 2 (4.3) 

i i 

where Hi and hi denote the SM and dark sector Higgs, respectively, and Qi and qi denote 
their SM and dark sector hypercharges. Integrating out the SM fields, Hi and Dy, generates 
a cross term eD y (Dy) in the low-energy theory. Thus, in the infrared, this induces an effective 

8 Strictly speaking, this is a form of gauge mediation according to the definition of Ref. [36] 
9 Kinetic mixing mediation is neglected in some recent (7(1) dark sector papers, for instance Ref. [37] focuses 
on a form of mediation that is sub-leading in e. 
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(4.1) 
(4.2) 



Fayet-Iliopolous D-term for D y , for which 

Vgaugc D eD y (D Y ) = tD y (4.4) 
£ = e (D Y ) = S^- cos 2/3 v 2 (4.5) 

where in the last equality we have substituted in for (Dy) from the MSSM. For e ~ 
1CU 3 — 10~ 4 , £ is at the GeV 2 scale. Thus, given the minimal assumption of kinetic mixing 
and SUSY, we obtain precisely the right scale to account for PAMELA and ATIC with the 
Sommerfeld enhancement. 

With an effective FI term at low energies, it is straightforward to break SUSY in the 
dark sector. In particular, a generic superpotential for the dark Higgses, hi, will break SUSY 
because the F and D terms cannot be simultaneously set to zero. While this SUSY breaking 
generates scalar soft masses, it does not generate soft masses for gauginos. Moreover, since 
SUSY is broken within the dark sector this will typically introduce a massless dark sector 
Goldstino 10 , assuming the absence of explicit SUSY breaking operators. So, within this 
mechanism there are light fermions. We present a concrete model of this type in Section 5, 
and mention the possibility of associated missing energy signals in Section 6.1.2. 

In the opposite extreme, we can take the superpotential to be less generic, or perhaps even 
trivial, and so the dark Higgs potential is dominated by D-terms. Here the dark Higgses simply 
align to set the dark hypercharge D-term to zero. In this limit SUSY is actually preserved 
in the dark sector, but a super-Higgs mechanism will generate a GeV scale dark sector that 
may still be consistent with a Sommerfeld enhancement and am z mass splittings for DM. A 
more minimal superpotential also can imply the existence of light pseudo-Goldstone fields and 
their superpartners. Finally, we note that unlike the SUSY breaking case, this super-Higgs 
scenario will also generate GeV scale gaugino masses. 

4.2 SUSY Breaking From 5 + 5 Messengers 

In this section, we elaborate on the gauge mediation proposal of [28] in which a multiplet of 5+ 
5 messengers is charged directly under both the dark sector and SM, thereby communicating 
SM SUSY breaking to the dark sector. We consider the additional possibility that the lightest 
component of the messenger supermultiplet is in fact the DM. 

Let us determine the various contributions which set the scale of masses for the scalar 
and fermion components of the DM 5 + 5 multiplet. First, we assume that the fermions have 

(2 3) 

a SUSY mass, ' , that splits the doublet and triplet components. Second, in the case 
where low-scale gauge mediation explains SUSY breaking in the MSSM, then the doublet and 
triplet scalars of the 5 + 5 receive identical soft mass contributions to that of the sleptons and 

(2 3) 

right-handed down squarks of the MSSM. We denote this contribution by mf' ' ~ 100 GeV. 
Finally, the scalar DM can in principle also receive soft mass contributions from whatever 
dynamics set its fi and B/j, terms, which we denote by B^ 2 ' 3 \ Instead of specifying these 

10 The gravitino will eat a linear combination of this field and the Goldstino associated with SUSY breaking 
in the MSSM 
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dynamics, we will choose a model-independent parameterization for the DM supermultiplet 
masses. The scalar doublets and triplets of the 5 + 5 have a scalar mass matrix given by: 



M 



2,3 



ill 



(2,3) 



+ 

5^(2,3) 



(2,3) 



m 



(2,3)' 



B/i (2,3) 
2 

+ 



(2,3) 



(4.6) 



whose eigenvalues, m±' 3 \ are given by 



in 



(2,3) 



+ m 



,(2,3) 



(2,3) 



(4.7) 



The lightest component of the doublet supermultiplet corresponds to dark matter, and we 
choose it to have mass ~ 500 — 800 GeV, which is favored by ATIC. 

Note that the messenger supertrace of the 5 + 5 mass matrix is non-zero, and proportional 

(2 3) 

to ms . As discussed in [38], this non-zero supertrace generates a logarithmically UV 
sensitive soft mass for the dark sector scalars. Indeed, since the messenger supertrace is 
positive, this implies a negative soft mass for the dark Higgs: 



f-V 



+ 3 



log 



3v 
m 2 f 



C a Sn 



(4.8) 



where C a is the dark scalar's quadratic Casimir, Sq is dark matter's Dynkin index, and Ajjv is 
set by the messenger scale of SUSY breaking to the SM. The negative soft mass squared allows 
for Gd ar k to break. This gives us a way to break the symmetries, independent of the effect of 
RGE running. It is our assumption that the contributions due to running are suppressed. For 
low-scale gauge mediation, Ajjv ~ 30 — 100 TeV, and because of this logarithmic enhancement 
and the combined effect of five messengers, we find that our desired scale of m\ ~ 1 GeV 2 
implies that m\ ' ' ~ 50 GeV. This indicates a bit of tension numerically because we expect 

(2 3) 

that ms ' is set by the SM soft mass scale of hundreds of GeV. 

Additionally, if we want fermionic DM, then there is the additional constraint that the 
fermion is the lightest component of the dark matter supermultiplet: thus (i?^ 2 )) 1 / 2 < 

(2) (2) 

m s <C rrij: . Since the DM B^i contribution breaks the dark sector R-symmetry, the gaug- 
ino soft masses are suppressed if we assume that the triplet component satisfies the same 
condition: 



-S 



2 7^ ~T~ 3~ 



rn 



(2) 



(3) 



(4.9) 

This implies light gauginos and the generic prediction is that fermionic dark matter implies 
that the lightest dark sector particle is a mostly-gaugino fermion. This conclusion can be 



(2) 
/ 



avoided by raising Bfi^ while maintaining (.B// 2 )) 1 / 2 <C m 

The dark sector Higgses require GeV scale fi and Bfi terms to help break dark gauge 
symmetry and lift runaway directions. These terms can be generated by additional dynamics 
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that communicate SM SUSY breaking to the dark Higgses, in general also resulting in new 
two-loop contributions to the dark scalar masses. We will assume that these contributions to 
m\ are subdominant to the usual gauge mediation contributions of Eq. 4.8. A recent paper 
identifies a class of general gauge mediation models that satisfy this assumption [39]. 

Let us note that while Eqs. 4.8 and 4.9 are approximations, our benchmark model of 
gauge mediation in Section 5.2.2 employs the full expressions of Ref. [38]. 

5. Benchmark Models 

In this section, we present four detailed benchmark dark sector models and their spectra. 
The models break dark gauge symmetry and custodial symmetry, generating the dark matter 
splittings and couplings necessary to explain the astrophysical data, as explained in Section 
2. These examples illustrate some of the theoretical issues discussed above, and their spectra 
and couplings serve as starting points for thinking about the types of cascades that can occur 
in GeV scale dark sectors. We begin in Section 5.1 with two non-SUSY models, where the 
GeV scale is put into the scalar potential by hand. We then consider two SUSY examples in 
Section 5.2, where the GeV scale is generated radiatively in the dark sector from interactions 
with the Standard Model. 

For each example we consider an SU(2) x U(l) y dark sector and triplet dark matter. We 
take the Majorana components of the dark matter fermions to be split by enough to avoid 
direct detection bounds, for example by one of the mechanisms discussed in Section 3. We 
then calculate the radiative splittings among the triplet, induced by dark symmetry breaking, 
as in Section 3.1. We take the ground state to correspond to dark matter, the heaviest excited 
state to allow for the XDM explanation of INTEGRAL, and the first excited state to allow 
for the iDM explanation of DAMA. We allow complex parameters to carry imaginary parts 
in order to avoid unbroken CP symmetry in the dark sector which may lead to stable states. 
This is not necessarily a problem and may actually have additional interesting signatures, 
but we'd like to keep the spectrum as general as possible for the present discussion. 

5.1 Non-SUSY Benchmark Models 
5.1.1 Non-SUSY 1: Two Doublets 

We begin with the two doublet model of Section 2.1, where h\ and I12 have dark quantum 
numbers 2_ 1 / 2 , 2i/ 2 - We have chosen a benchmark which breaks charge and radiatively 
generates the XDM and iDM splittings, and we have calculated its mass spectrum (Fig. 5). 
As discussed in section 2, the custodial SU(2) symmetry of the Higgs sector determines the 
tree-level gauge boson spectrum. The gauge bosons that couple between the different dark 
matter mass eigenstates, w±, do not mix with the b and are degenerate in mass. Custodial 
symmetry is broken at one-loop and in general due to higher-dimensional operators. The 
DAMA inelastic scattering is therefore suppressed relative to models where custodial sym- 
metry is broken at tree-level. For this benchmark, we induce the iDM coupling by including 
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the dimension 6 custodial-breaking operators cJ\hiDhi\ 2 and | /i2^^2 1 2 ? with coefficients 
cj and C2 expressing the loop-suppression. 

For the benchmark, we choose the gauge couplings: g = 0.46 and g y = 0.19. The dark 
matter hypercharge is chosen to be y^m = 1/2. In the limit of small charge breaking, this 
choice leads to one small and one large dark matter splitting, as discussed in Section 3.1. In 
terms of the potential of Eq. 2.3, the parameters are: v± = 1.5 GeV, V2 = (1.5 + 3.2 i) GeV, 
Ai = 0.5, A2 = 0.3, A3 = —0.031, A4 = 0.5, and cos a = 0.8. The coefficients of the 
custodial-breaking dimension 6 operators are chosen to be cj = 2.8 x 10" 4 GeV~ 2 and 
c t = _ 5 .7 x iQ-4 GeV" 2 . 
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Figure 5: The spectrum of Non-SUSY 1, our two doublet non-SUSY benchmark. The left side 
shows the radiative mass splittings of the components of the dark matter triplet, measured from the 
ground state. The splittings allow for the XDM and iDM explanations of INTEGRAL and DAMA, 
respectively. The right side displays the spectrum of the GeV-scale dark sector. The b fractions of 
the gauge bosons are indicated and determine how strongly each gauge boson couples to Standard 
Model electromagnetic current. Because of custodial SU(2), two of the gauge bosons are degenerate 
and do not mix with the b at tree-level, and these are the gauge bosons that couple between different 
dark matter states. They do mix with the b at one-loop, inducing a suppressed iDM coupling, and 
we include the dimension 6 operators c\\h\Dh\\ 2 and c^\hiDh-i^ in order to parametrize custodial 
breaking corrections. The parameters of this benchmark are listed in the text. 



5.1.2 Non-SUSY 2: Two Doublets and One Complex Triplet 

We now add a complex triplet Higgs <1> to the two doublet model, with dark quantum numbers 
3q. The triplet vev breaks custodial symmetry, causing all gauge bosons to mix with the b 
and inducing the iDM coupling at tree-level. We take the triplet to be complex. While not 
the minimal possible choice, it has a more straightforward SUSY extension. We again choose 
a benchmark that breaks charge in the doublet sector and radiatively generates the XDM 
and iDM dark matter splittings. We have calculated its spectrum (Fig. 6). 
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For this benchmark, we choose the gauge couplings g = 0.23 and g y = 0.75, and the dark 
matter hypercharge is chosen to be ydm = 0.3. The potential is similar to Eq. 2.6 except we 
take to be complex: 



field 


\ hi | h 2 




charge 


2-1/2 | 2]y 2 


3o 



V(h u hz,*) 



V(h 1 ,h 2 ) + ^-(\r 



$t$ 



(5.1) 
(5.2) 



cih{®hi + c 2 h\§h 2 + c 3 hi e®h 2 + h. 



where the first term is the two doublet potential of Eq. 2.3. For the two doublet sector 
we choose the parameters: v\ = 1.8 GeV, v 2 = (1.8 + 1.4 i) GeV, Ai = 0.71, A2 = 0.47, 
A3 = 0.33, A4 = 0.099, and cos a = 0.052. The parameters involving the triplet are chosen to 
be: v$ = (1.1 + 0.61*) GeV, A$ = 0.51, c x = (0.054 + 0.47 %) GeV, c 2 = (0.74 + 0.69 i) GeV, 
and c 3 = (0.61 + 0.81 i) GeV. 
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Figure 6: The spectrum of Non-SUSY 2, our two doublet and one complex triplet non-SUSY 
benchmark. The left side displays the radiative mass splittings of the dark matter triplet, measured 
from the ground state. The splittings can account for the XDM and iDM explanations of INTEGRAL 
and DAMA, respectively. The right side shows the dark sector spectrum, and the b fractions of the 
gauge bosons are indicated. The triplet vev breaks custodial SU(2), and all 4 gauge boson mass 
eigenstates mix with the b at tree-level, although for this example most of the b is contained in two of 
the mass eigenstates. The parameters of this benchmark are listed in the text. 



5.2 SUSY Benchmark Models 

Now we consider two SUSY benchmarks, where the GeV scale is generated radiatively from 
interactions with the Standard Model, as discussed in Section 4. Since our models all employ a 
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kinetic mixing, they all receive SUSY breaking contributions from kinetic mixing mediation, 
as discussed in Section 4.1. In our first example, SUSY 1, this the only source of SUSY 
breaking, however as noted in appendix B, it is difficult within this framework to break 
charge with only one hypercharge neutral triplet. We circumvent this in this example by 
adding a second complex triplet and taking the triplets to have dark hypercharge. For our 
second example, SUSY 2, we add an additional gauge mediation source for GeV-scale SUSY 
breaking by taking dark matter to be charged as a 5 + 5 of the SM. Dark matter then acts 
as a messenger of gauge mediation, as discussed in section 4.2. For this setup, we can break 
charge with two doublets and one hypercharge neutral triplet. 

5.2.1 SUSY 1: Kinetic Mixing Mediation with Two Doublets and Two Triplets 

For this benchmark, we have two doublets, hi and h 2 , and two complex triplets, <&i and & 2 , 
with dark quantum numbers 2_ 1 / 2 , 2^2, 3i, and 3_i. We have chosen triplet hypercharge as- 
signments that allow Yukawa couplings between doublets and triplets, otherwise there may be 
pseudo-Goldstone bosons in the spectrum. The GeV scale is generated in the dark sector from 
kinetic mixing, as described in Section 4.1. The most general renormalizable superpotential 
for two doublets and two triplets with these charge assignments is: 



field 1 hi 


h 2 | 


$1 [ 


$2 


charge || 2_i/ 2 


2l/2 


3i 


3-i 



W = Hhh\eh 2 + /U$Tr [$i$ 2 ] + Xyh^eQxhi + X 2 h^e<^ 2 h 2 (5.4) 

We include GeV scale fx and B/j, terms for the doublets and triplets because they help break 
the dark gauge symmetry and lift runaway directions. We do not include the effects of 
running from the TeV scale to the GeV scale, which we take to be subdominant. Kinetic 
mixing mediation already gives negative scalar soft mass squareds at tree-level, leading to the 
breaking of dark gauge symmetry. 

We include a kinetic mixing coefficient of e = 2 x 1(U 4 , in terms of Eq. 2.1, which auto- 
matically generates the GeV scale in the dark sector. Our benchmark radiatively generates 
the XDM splitting, but unfortunately the smaller dark matter splitting is too large to account 
for iDM. We have calculated the mass spectrum (Fig. 7). The triplet vevs break custodial 
SU(2) at tree-level and all gauge bosons mix with the b. The gauginos and Higgsinos are 
strongly mixed after dark symmetry breaking, but for this example the lightest fermion is 
a mostly gaugino-like Goldstino with a mass of only ~ 2 MeV. Such a field is present be- 
cause SUSY is broken within the dark sector itself. The second lightest fermion, with mass 
~ 190 MeV, is lighter than the lightest gauge boson. Thus, the dark gauge bosons will cascade 
into these light fermions, rather than SM lepton pairs. The second lightest fermion decays to 
the lightest fermion and a SM lepton pair through a 3-body decay, which can account for the 
astrophysical lepton production and lead to visibly displaced vertices at colliders. Another 
possibility, not realized in this example, is to have an approximately supersymmetric dark 
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sector, with a kinetic mixing mediation induced super-Higgs mechanism at a GeV. Gauginos 
then reside in massive vector supermultiplets and get GeV scale masses. 

For this benchmark, we have chosen the gauge couplings g = 0.22 and g y = 1.2, and 
dark matter hypercharge ?/dm = 1/5. The superpotential Yukawa couplings are chosen to be 
Ai = 1.7 + 0.022 % and A 2 = 0.5 + 1.8*. For the doublets we choose n h = (0.11 + 0.63 *) GeV 
and (Bfj) h = (0.74 + 0.69*) GeV 2 . For the triplets we choose /x$ = (0.51 + 0.83 i) GeV and 
(Bp)* = (0.57 + 0.59 i) GeV 2 . 
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Figure 7: The spectrum of SUSY 1, our two doublet and two complex triplet SUSY benchmark 
with kinetic mixing mediation. The left side displays the radiative mass splittings of the dark matter 
triplet, measured from the ground state. The larger splitting allows for the XDM explanation of 
INTEGRAL, but the smaller splitting is too large to explain DAMA with iDM. The right side shows 
the dark sector spectrum, and the b fractions of the gauge bosons are indicated. The triplet vevs 
break custodial symmetry, and all four gauge boson mass eigenstates are part b at tree-level. The 
gauginos and Higgsinos are strongly mixed after dark symmetry breaking. The lightest fermion, with 
mass ~ 2 MeV, is mostly gaugino and light because gauginos get no soft masses from kinetic mixing 
mediation. The second lightest fermion has a 3-body decay to the lightest fermion and a SM lepton 
pair, which can account for astrophysical lepton production and lead to a visibly displaced vertex at 
colliders. 



5.2.2 SUSY 2: Gauge Mediation with Two Doublets and One Triplet 

For this benchmark, we supersymmetrize the Higgs content of our Non-SUSY 2 benchmark, 
including SUSY breaking contributions from both kinetic mixing mediation and gauge medi- 
ation with dark matter messengers. The most general renormalizable superpotential for two 
doublets and one triplet is the following: 
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As in the SUSY 1 benchmark, we include GeV scale fi terms for the doublets and triplet 
and do not include the effects of running from the TeV scale to the GeV scale. There are 
already negative scalar soft mass squareds at tree-level because of the nonzero dark matter 
supertrace [38], leading to the breaking of dark gauge symmetry. For this example, it is not 
necessary to include GeV scale B[i terms for the doublets or triplet. 

We have chosen a benchmark which generates a GeV scale dark sector with charge break- 
ing and custodial breaking, and which leads to radiative XDM and iDM splittings. We have 
calculated the mass spectrum (Fig. 8). The gauginos and Higgsinos are strongly mixed after 
dark symmetry breaking, but the three heaviest fermions with masses near ~ 5.5 GeV are 
almost pure Higgsino mixtures. The spectrum is slightly split by a small separation between 
the dark jjl and soft mass scales. The gauge couplings are chosen to be g = 0.3 and g y = 0.37 
and the dark matter hypercharge is ydm = 1/2. The kinetic mixing is e = 7 x 10 -5 in terms 
of equation 2.1. The messenger scale of SUSY breaking to the Standard Model, which enters 
the log divergence of Eq. 4.8, is chosen to be Ajjv = 30 TeV, corresponding to low-scale gauge 
mediation. The standard model doublet dark matter mass components, in terms of Eq. 4.6, 
are given by mf* = 800 GeV, = 50 GeV, and Bfi^ = (40 GeV) 2 . As discussed in 

section 4.2, the small soft mass is needed to generate the GeV scale in the dark sector, and 
we choose a small Bfi, keeping dark matter fermionic. For the standard model triplet compo- 
nents of the dark matter 5 + 5, we choose the parameters nij = 840 GeV, m s = 50 GeV, 
and Bfi^ = (300 GeV) 2 . The larger B[i for the triplet leads to GeV scale gaugino soft 
masses in the dark sector (see Eq. 4.9 and the surrounding discussion). The superpotential 
parameters are fj, h = (0.27 + 0.28 i) GeV, /i$ = (2.52 + 3.48 i) GeV, and A = 0.29 + 1.51 i. 

6. Signals of a Non-Abelian Dark Sector at the Tevatron and the LHC 

In this section we discuss the collider phenomenology associated with the models presented 
in the previous sections. In the first part of this section we analyze the generic predictions 
associated with a non-abelian dark sector that is linked to the SM only via kinetic mixing. In 
the second part we present the signals expected in supersymmetric versions of such models. 
Throughout, we limit the discussion to the Tevatron and LHC. It is important to realize that 
in the case of a GeV scale dark sector such high-energy accelerators are needed more for their 
luminosity than their energy reach, which is considerably higher than the dark sector scale. 
In supersymmetric implementations, colored MSSM superpartners can be copiously produced 
at hadron colliders. Their subsequent decay into dark states can produce spectacular signals 
involving multiple lepton jets. We leave it for future work to investigate the phenomenology of 
these models at low-energy experiments, but sec Ref. [11] for low-energy signatures of similar 
models. 

A new sector of light particles with very weak couplings to the Standard Model have 
been discussed in detail in the context of the "Hidden Valley" models [40]. Their collider 
phenomenology was investigated in [41, 42]. In particular, the modifications such models can 
introduce to the decay chains of the MSSM was clarified in Ref. [43]. Here, we focus on the 
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Figure 8: The spectrum of SUSY 2, our two doublet and one complex triplet SUSY benchmark with 
both kinetic mixing mediation and dark matter messengers charged under the Standard Model. The 
left side displays the radiative mass splittings of the dark matter triplet, measured from the ground 
state. The splittings can allow for the XDM and iDM explanations of INTEGRAL and DAMA, 
respectively. The right side shows the dark sector spectrum, and the b fractions of the gauge bosons 
are indicated. The triplet vev breaks custodial SU(2) and three of the gauge boson mass eigenstates 
are part b at tree-level. The dark spectrum now includes GeV scale fermions, and the gauginos 
and Higgsinos are strongly mixed after dark symmetry breaking. The three heaviest fermions, with 
masses near ~ 5.5 GeV, are almost pure Higgsino mixtures. The spectrum is slightly split by a small 
separation between the dark [i and soft mass scales. The parameters of this benchmark are listed in 
the text. 

particular scenario which uses the kinetic mixing as the essential link between the SM and the 
dark sector. In addition, motivated by astrophysical observations, we allow the dark sector 
to decay back to light leptons (e ± and /r* 1 ) only. For the purpose of this paper, we do not 
concern ourselves with a possibly small branching fraction into pions. 

The pair production of the dark matter states at colliders is certainly possible if they 
happen to carry SM weak charge. However, their detection proves extremely difficult since 
they are not accompanied by any hard object. Even if the excited states of its SM multiplet 
are produced, their decays are too soft to trigger on since they are separated by only ~ GeV 
(notice that this splitting is generated by the SM gauge interactions and are of order aMz 
[34]). 

6.1 Production and decay of dark gauge bosons and Higgses 

As discussed in detail in appendix A, the kinetic mixing induces two important, e suppressed, 
couplings: The SM electromagnetic current is now also charged under the dark gauge bosons; 
the SM Z° boson is now coupled to the dark hypercharge current. Before discussing each of 
these couplings and their impact on collider signals, let us briefly discuss the decays of the 
dark gauge boson and the dark Higgses. 
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6.1.1 Dark gauge boson and Higgs decay chains 




Figure 9: Typical decay chains starting with a dark gauge boson, labelled 7' in this plot. The dark 
decay chain can have several stages and involve additional dark sector states, such as other dark gauge 
bosons (labeled w',z' in this figure), and dark Higgses (labelled h'). 

The non-abelian nature of the dark sector implies the presence of complicated decay 
chains. Some of the typical decays chains are shown in Fig. 9. In the dark sector, gauge boson 
mass eigenstates are generically mixtures of all four SU(2)daik x U(l) y gauge eigenstates. In 
Fig. 9 and the rest of this section, we have used 7' (and also w' and z' in this figure) to 
denote any one of these mass eigenstates. For an abelian dark sector with kinetic mixing 
with the SM, 7' decay leads to a di-lepton final state, shown in the first panel from the left 
of Fig. 9. On the other hand, a non-abelian dark sector, like one of the examples considered 
in this paper, leads to complicated decay chains, such as the ones shown in the rest of Fig. 9. 
The dark Higgs sector, necessary to break the non-abelian group, may also participate in 
such cascades as shown in the right two panels of Fig. 9. Such cascades inevitably produce 
multiple, easily > 2 and possibly 8, final state leptons, which provides a unique signature of 
the non-abelian nature of the dark sector 11 . We expect the decay between dark states to be 
generically prompt. Therefore, the decay length is dominated by the very last decays back 
into SM leptons. A rough estimate for a generic decay is then, 

^J^-,r,ro- -(^)(^)". 

With moderate boost 7 ~ O(10), this may lead to a displaced vertex if e < 10 -4 . 

To be observable at hadron colliders, the dark boson which initiates such a cascade must 
carry px ~ O(10s) GeV. Therefore, regardless of the precise nature of the cascade which 
ensues, its decay products have small opening angles 59 ~ my/pr < 0.1. Those decay 
products will eventually decay into several collimated SM leptons. A collection of more than 
2 hard and collimated leptons is dubbed a "lepton jet" [28]. 

6.1.2 Displaced vertices and missing energy 

While Eq. (6.1) is the generic estimate for the resulting decay length of dark cascades, 
there are several exceptions which may result in more noticeably displaced vertices or missing 
energy in lepton jets. 

11 Sometimes phase space constrains the flavor of the lepton. For example, a GeV dark gauge boson cannot 
decay into more than 4 muons 
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If it is kinematically forbidden for a dark gauge boson to have 2-body on-shell decays 
within the dark sector, then the dark gauge boson may decay directly into two leptons. 
However, a noticeable exception occurs when the 3-body decay 7' —* a'*b'i — > b^b^b'^ is kine- 
matically allowed, where a', b' can be either dark gauge boson or dark Higgs states. In this 
case, there is an additional suppression of (5m/m a i) 5 x (3-body phase space) on the decay 
width, where 5m ~ my — £^ my , and we have used my ~ m a i in this estimate. This decay 
channel can be competitive and even dominate over the direct decay into 2 leptons. In par- 
ticular, when the decay into SM leptons is strongly suppressed (dark pseudoscalar decay) or 
all together forbidden (dark fermion decay), the 3-body process may dominate and lead to a 
displaced vertex. The impact parameters of multiple leptons associated with this displaced 
vertex will not be correlated with each other since they come from the decays of different 
resonances b' 12 3 . 

If the lightest dark sector state is a dark Higgs, h' , it cannot directly decay into SM 
leptons (unless it mixes the SM Higgs, see Ref. [19, 1]). In this case, the dark Higgs will 
either decay into 4 leptons through two off-shell gauge bosons, shown in the left panel of 
Fig. 10, or into 2 leptons through a one-loop decay. Either way, such a decay leads to a 
very long life-time, cr ~ 0(km) for m^ < GeV. In this case, dark cascades which involve 
this lightest scalar contain missing energy as it escapes the detector. These cascades can still 
produce observable lepton jets because, in addition to missing energy, one still gets leptons 
from the intermediate steps of the decay, such as h! i — > a'h' followed by a' — ► lepton pairs. 
In this case, the lepton jet contains missing energy that is collimated with the leptons of the 
same cascade. 




Figure 10: Two possible decay channels if the lightest dark sector state is a scalar from the dark 
Higgs sector. 

An additional source of missing energy comes in a supersymmetric dark sector with Im- 
parity. The lightest dark supersymmetric particle (LDSP) may be stable if the gravitino 
is heavier. Otherwise, it may eventually decay into the gravitino. Either way, it carries 
with it missing energy. Unless the MSSM sector decays directly into the LDSP, in which 
case there may be no lepton jets, missing energy due to the LDSP will be collimated with 
the visible lepton jets, very similarly to the non-SUSY case. Such correlations provide an 
additional handle on the reconstruction of these events since we know the direction of the 
missing particles and can treat them as having vanishing masses. We provide an example of 
such a reconstruction in the case of rare Z° decays below. 
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6.1.3 Direct Production 




Figure 11: Direct production of a dark gauge boson in a process very similar to prompt photon 
production in the Standard Model. 




Figure 12: In the left pane, we show the rate of direct production of the dark gauge boson as a 
function of e e g/e, where e c ff is the effective coupling of dark gauge boson to the Standard Model 
fields. 



The kinetic mixing between the dark force carrier and the SM photon induces a small 
dark charge for electromagnetically charged SM fields. Consequently, the dark gauge boson 
can be directly produced in colliders via a process analogous to prompt photon production 
in the SM, shown in Fig. 11. 

In the left panel of Fig. 12, we present the production rate of dark gauge bosons as a 
function of e c g/e, where e e g = eecos 9wfb is their effective gauge coupling to SM fields 12 and 
fb is the fraction of the dark hypercharge gauge boson in a given dark gauge boson mass 
eigenstate. In the right panel of Fig. 12 we plot the inclusive differential cross-section of dark 
photon (7') production at the LHC and the Tevatron with e e g = 10 _3 e. 

12 The simulation was actually of prompt photon production with PYTHIA [44] and the resulting cross- 
section was multiplied by a factor of e^/e 2 . 
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Figure 13: pt distributions for cascades resulting in 4 (left) and 8 (right) leptons, for events with 
p T > 50 GeV for 7'. 
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Figure 14: The differential cross-section as a function of the pr of 7' at the LHC (y/s = 14 TeV) 
after including muon triggers, demanding either a single muon with px > 7 GcV or two muons with 
Pt > 3 GeV. Proper r\ cuts were imposed and each event was required to contain at least 3 leptons. 



After dark vector bosons are produced, they typically cascade down to multiple leptons 
that form a lepton jet as discussed above. These leptons carry a significant amount of pr, as 
shown in Fig. 13. At CMS, the Level 1 Dimuon trigger (2 muons with pt > 3 GeV in \r/\ < 2.4) 
or single muon trigger (1 muon with pj- > 7 GeV and \rj\ < 2.4) should be able to detect those 
events that contain muons [45]. The electron triggers are single e (isolated Et > 26 GeV), 
double e (isolated Et > 14.5) and double relaxed e (not isolated Et > 21.8 GeV). Since the 
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Figure 15: The fraction of events with a maximum opening angle 9 between leptons in a lepton jet, 
which contains 4 leptons. 

resulting electrons are unlikely to be isolated "electromagnetic" objects, the double relaxed 
e is probably necessary. We will conservatively assume that muons alone are triggered on. 
In Fig. 14, we show the differential cross section of dark 7', taking into account the simple 
requirements on muon triggering. 

6.1.4 Distinguishing Leptons 

Let us discuss the issue of discriminating individual leptons within a given lepton jet 13 . In 
our present discussion, we focus on muons. In Fig. 15 we plot the maximal opening angle 
between any two of the four leptons. At such high momenta, the resulting decay products are 
highly collimated with an initial opening angle of approximately 9 ~ my/pT < 0.1, which 
can be as small as 10~ 2 . By the time these muons reach the first layer of the muon system, 
they typically acquire a sufficient separation to be distinguished. For example, as depicted in 
Fig. 13, a typical scenario will have two muons with average pr ~ 20 GeV and Apx ~ 5 GeV. 
Without even including the initial lepton jet opening angle, we estimate that the acquired 
separation is about 10 cm (in the CMS detector), which is greater than the cell size of ~ 4.5 
cm. The separation between two same sign muons is proportional to Apx / '(pj. 1 Pj? ') ■ For 
a given lepton jet pt, since both Apt and pt are inversely proportional to the number of 
leptons, higher multiplicities actually result in larger separations. We also notice from Fig. 13 
that leptons in lepton jets typically have different pt, such that Apx/pT ~ 20% or so. LHC 
detectors can achieve better muon momentum resolution. For example, CMS can achieve 
Apr/pr ~ 1% (about 10% with muon system only) in the momentum regime of interest 

13 We are grateful for valuable discussions with Jim Olsen on this subject. 
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[46]. ATLAS can achieve a similar precision [47]. Finally, the muon isolation separation 
defined by CMS can be as small as AR = 0.01. The angular resolution is even smaller, 
about 2 mrad[46] 14 . Thus, it is reasonable to assume that CMS will be capable of resolving 
several, if not all of the muons. The primary background arises from K and tt decays and 
J/ip — * (muons coming off soft jets can be vetoed with isolation cuts), and possibly 

from other heavy flavor decays. The high lepton multiplicity in those events and the lack of 
hadronic activity around the lepton jet should be sufficient to fight these backgrounds and 
obtain a clean sample. However, a more careful collider analysis is certainly warranted, but 
is beyond the scope of the present work. 

6.1.5 Rare Z° decay 

As discussed in Appendix A, the kinetic mixing also induces a coupling eZ^Jjf, where is 
the dark hypercharge current. Thus, we can produce dark hypercharged states through rare 
decays of the Z°, shown schematically in Fig. 16. The e 2 suppression makes LEP searches 
irrelevant due to luminosity limits, but the Tevatron and LHC may probe such events. The 
decay branching ratio to any particular dark sector state d{ can be written as 



BR(Z C 
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where depends on decay matrix element and is proportional to the number of degrees 
of freedom of d^. The total branching ratio into the dark sector will scale linearly with the 
number of dark sector states, which could be easily O(10) in our scenario. 
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Figure 16: Z° production and two possible decay channels into the dark sector. On the left we 
depict a decay into the dark Higgs sector. Fermionic channels, such as the one shown on the right, 
arc dominantly associated with the Higgsino states possible in supersymmetric versions of the model. 



The SM photon does not couple to the dark sector states. However, there is a "continuum" 
contribution to the same amplitude through off-shell dark photon, qq — > 7'* — > d^di, which 
is proportional to e 2 ff oc e 2 . Depending on the spectrum and couplings in the dark sector, it 
could have important contributions to the signal off the Z-peak. In this section, we will focus 
on the contribution within the Z resonance. 

14 This is the resolution quoted for a single hit. It is beyond our abilities to evaluate how the resolution 
deteriorates with multiple collimated muons. 
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The production rates of dark sector states at the Tevatron and LHC are shown in Fig. 17 
[48]. We present rates coming from decays into bosonic (denoted by h') and fermionic (denoted 
by /') dark sector states. In the context of the SUSY models discussed later in this section, 
these bosonic and fermionic states could refer to dark Higgs bosons or Higgsinos, respectively. 
On the other hand, the collider phenomenology is similar if other possible dark sector states 
decay into lepton jets. A cut of \rj\ < 2.4 has been imposed on the direction of the lepton 
jets. The difference in rates between the fermionic and bosonic channels results from the r\ 
cut, the boost of the Z° in the lab frame, and the fact that fermions are more likely to be 
emitted along the boost direction because of angular momentum conservation. 




Figure 17: Left: The production rate as a function of the branching ratio of the decay: Z° — ► dark 
states. The solid and dashed lines are for Z° decays into dark sector scalars and fermions, respectively. 
Right: lepton jet px distribution resulting from Z° decays. 

As can be seen from Fig. 17, the lepton jets produced in this way are peaked towards 
^lepton jot ^ 0.5 M^. Therefore, they are typically harder than the lepton jets resulting from 
the prompt production of dark gauge bosons. As we have discussed in Section 6.1.3, such 
harder lepton jets will be easier to trigger on. However, we expect the efficiency of identifying 
different leptons in a lepton jet will be lowered as it is oc 1/pt- 

Reconstructing the Z° is not difficult and helps to reduce the background. With enough 
statistics, it is even possible to study the angular distribution of the resulting lepton jets and 
get a handle on the spin of the dark sector states as demonstrated in Fig. 18. About 5000 
events are used in this plot. We see that the expected rise for cos 9 ~ ±1 from the fermionic 
decay channels is washed out due to the large boost of Z° and the |r/| cut. However, the 
resulting distribution is still quite different from that of the bosonic decay channel. 

6.2 Collider signals of supersymmetric models 

In this section, we discuss the collider signals associated with supersymmetric models. In 
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Section 4, we have focused on models with low supersymmetry breaking scale. However, the 
present discussion of the resulting collider signatures is largely independent of that scale or 
other MSSM details since we will not consider any specific superpartner spectrum. In that 
sense, models with higher supersymmetry breaking scales, such as the Planck slop option 
suggested in Ref. [28], are only different in that the gravitino is heavier. Hence, the end of the 
dark sector decay chain will not involve the gravitino. However, this does not have a visible 
effect on the collider signals. Even in the low scale models where the gravitino is light, the 
decay length of the dark sector LSP to the gravitino is much larger than the detector size. 
That said, it is important to note that the collider signatures discussed in this section are 
based on the assumption that the MSSM LSP dominantly decays into the dark sector. 

With supersymmetry, the dark sector states are dominantly produced from cascade decays 
of MSSM colored superpartners, such as gluinos and squarks. These particles follow typical 
MSSM decay chains down to the MSSM LSP (not the gravitino). The effect of the GeV dark 
sector is to extend and/or modify the decay chains following MSSM LSP production [28]. We 
begin by summarizing the main features of such cascades. 

Let us first note, however, that a notable exception occurs if dark matter is part of a 
pair of 5 + 5 under the SM gauge groups. An example of such a model was presented in the 
benchmark of Section 5.2.2. The rate for the production of the colored components of such a 
pair is shown in Fig. 19. Thus, the LHC has great potential for producing such states up to 
about 2 TeV. As pointed out in Ref. [28], as long as the colored particles decay only through 
higher dimensional operators they will be long-lived and may have decays with very distinct 
signatures [49]. We will not elaborate on these possibilities but refer the interested reader to 
the detailed discussion in Ref. [28]. 
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Figure 19: Production rate of one set of 3 + 3 vector-like quarks, which can be part of the dark 
matter multiplet. The rate should be scaled by the number of such representations and the dimension 
of the dark matter representation under Gdark- 

6.2.1 MSSM decays into the dark sector 

Kinetic mixing implies that if the MSSM LSP is a neutralino then it decays into dark sector 
states with a lifetime of 

T LSP-vh+h ~ ( a f rl 7le 2 ^LSp) 

where fg is the bino fraction of the MSSM LSP. In the low-scale gauge mediation models 
constructed earlier in this paper, it is possible for the gravitino to be significantly lighter than 
the MSSM LSP. When the gravitino is light, another possible decay channel for the MSSM 
LSP is LSP— > XsmG, where Xsm can be a photon, Z, or Higgs, depending on the model 
parameters and phase space. The decay lifetime can be estimated as 

/M I 5 SP \" 1 „ /100GeV\ 5 / y/F V 

r^,z M a ~ ) = 3.3 X 10- . (-g— ) J . (6.4) 

We see that the LSP dominantly decays into the dark sector instead of the gravitino. However, 
the two channels can be comparable in certain regions of parameter space, such as ~ 0.1 
and a low supersymmetry breaking scale yF ~ ^Jrn^Mp ~ 10 TeV. 

When the MSSM LSP is a sfermion (£ or q), things become more subtle. One possible 
decay channel is through an off-shell gaugino with a significant bino fraction, / — > / + \* ~^ 
f + [dark sector states], shown the left panel of Fig. 20. Its decay lifetime can be estimated 
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Figure 20: Decay of sfermion LSP. 
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where Cf x is the effective fermion-x coupling, and P(mj/M x ) is a function with the limit 
P — > (rrir/M^ 4 for M x 3> m?. Another possible decay channel is / — ► / + b, shown in 



the right panel of Fig. 20, where b is the dark bino. However, as explained in Appendix A, 
in addition to kinetic mixing, this coupling has an additional suppression of order Mr/Mg. 
Hence, its lifetime is 
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Notice that when the off-shell gaugino state is dominantly bino, we have 
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Therefore, these two channels can be either quite different or comparable, depending very 
sensitively on the details of the model. In principle, these two decay channels are distinguish- 
able experimentally, as the three (two) body decay gives rise to three (two) different lepton 
jets, respectively. Notice also that in this case, it is easier for the channel that decays into 
the gravitino to be competitive as well, if F is close to tens of TeV. 

Once the cascade has progressed into the dark sector, it will decay through the mass 
hierarchy there. Several typical cascades in the supersymmetric dark sector are shown in 
Fig. 21. For example, the center panel on the second row is important when the fermionic 
dark superpartners are lighter than the dark gauge boson. An example of such a scenario is 
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Figure 21: Typical SUSY dark sector decay chains. The dark sector states, 7', h! and z' , can be 
either on or off shell. They will cascade further to produce lepton jets, similarly to the non-SUSY 
case. 



discussed in Section 5.2.1. All of the decay products in the same chain are collimated into 
one lepton jet with typical pt ~ 100 GeV. Therefore, it is difficult to uncover the details of 
the decay chain that produces a given lepton jet. 

Finally, we discuss the endpoint of the dark sector decay chain. First, consider the 
situation where the lightest dark sector particle is the LDSP, as in the benchmarks of sections 
5.2.1 and 5.2.1. If the gravitino is heavier than the LDSP, then the decay chain will end there 
with the LDSP escaping the detector, producing missing energy along the direction of the 
lepton jet. If the gravitino is lighter than the LDSP, as in the case of low scale supersymmetry 
breaking models, the last step of the cascade will be LDSP^ X^mG, where Xq,m corresponds 
to a light SM particle, such as a photon or lepton. For simplicity, we consider the case where 
the LDSB is mostly dark bino. We can estimate its decay lifetime as 



e 2 M- 5 
167T F 2 "" 



3.3 x 10 J s 




(6.8) 



which is clearly not relevant on the collider timescale. Second, we consider the case where 
there is also a dark sector gauge boson, b, that is lighter than the LDSP. This situation 
is not realized in our benchmarks of Section 5.2, but is certainly a possibility. In this case, 
the LDSP can decay through the channel, b — > bG, where b subsequently decays to leptons. 
There is no e 2 suppression here, but setting e to 1 in Eq. (6.8) gives a decay length inside 
the detector only for y/F < 10 TeV. Therefore, we can effectively think of the LDSP as the 
endpoint of the dark sector decay chain for most of parameter space. 

6.2.2 Extended discovery reach for direct electroweak-ino production 

The direct production of electroweak-inos is an important channel since it is independent 
of the existence of colored superpartners and may provide additional information on the 
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Figure 22: Pair-production of the SM LSP can result in spectacular lepton jet + MET events. On 
the left we depict the event topology. On the right we show a schematic representation of the resulting 
geometry. 




Figure 23: Left: The cross-sections for electroweak-ino production at the Tevatron. We have included 
both LSP pair production and, in the case of wino and Higgsino LSP, the production of closely 
degenerate states, as a function of M x . We choose the squark mass to be 750 GeV. Right: the 
fraction of events with 3 and 4 lepton jets within the central region |ry| < 2.4. 



properties of those electroweak-inos. In the conventional MSSM, it is usually difficult to 
see events with direct pair-production of electroweak-inos. In the case of direct MSSM LSP 
production, one has to trigger on some additional hard radiation, which has a lower rate 
and a large background. The pair-production of heavier electroweak-ino states which cascade 
down to the LSP may be easier to observe but suffers from large SM background. However, 
in the scenario we consider, the LSP of the MSSM, which we denote by Xo> wm decay further 
into the dark states [43] whose decays result in leptons and missing energy [28]. Such events 
are easy to trigger on since all the leptons carry significant amounts of pr- Since Xo is 
produced almost on threshold, its boost factor is order unity and the opening angle in the 
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Figure 24: Left: The cross-sections for electroweak-ino production at the LHC. We have included 
both LSP pair production and, in the case of wino and Higgsino LSP, the production of closely 
degenerate states, as a function of M x . We choose the squark mass to be 750 GeV. Right: the 
fraction of events with 1, 2, 3 and 4 lepton jets within the central region \r)\ < 2.4. 



decay Xo —> h DM XDM is fairly large. The resulting event geometry is striking and is depicted 
schematically in Fig. 22). 

In the left panel of Fig. 23, we show the rate of electroweak-ino pair production at the 
Tevatron. In the case of pure wino-like and Higgsino-like LSP, we have also included the 
production of the closely degenerate charginos and neutralinos. We then decay each LSP 
into a pair of lepton jets and study their kinematics. At the Tevatron, the neutralinos and 
charginos produced from qq' initial states are expected to have small boosts. Therefore, the 
majority of the resulting lepton jets are expected to be very central as shown in the right 
panel of Fig. 23, where we have required \rj\ < 2.4 for the lepton jets. In addition we see that 
the majority > 90% of the events have 4 lepton jets within the central region as illustrated in 
the left panel of Fig. 23. Since the presence of such lepton jets greatly enhances the possibility 
of triggering on such events and separating them from the background, we estimate a reach 
of about 300 GeV for pure Higgsino or wino LSP at the Tevatron. The case of pure bino is 
still difficult because of the suppression in rate. 

We have shown a similar study for the LHC in Fig. 24. At the LHC, the qq' initial state 
will carry significant boost. Therefore, as can be seen in the right panel of Fig. 24, there is a 
significant fraction of events with 3 or less lepton jets in the central region, especially for the 
smaller electroweak-ino mass M x < 400 GeV. On the other hand, as M x increases, the effect 
of the boost quickly decreases and the fraction of events with 4 lepton jets increases. Such 
lepton jets will give the LHC the amazing ability to probe bino production up to Mg ~ 1 
TeV, and wino or Higgsino production up to 2 TeV. 

6.2.3 Measuring the mass of the MSSM LSP 

In the sorts of SUSY events shown in Fig. 22, it is possible to use lepton jets for a measurement 
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Figure 25: Forming the invariant mass of all the lepton jets in the events can lead to a measurement 
of the mass of xo> shown here as an edge at M Xo . Incorrect pairings of the lepton jets are included. 
However, we assumed that individual leptons are properly bunched with the correct lepton jet. 

of the mass of xo- There are two lepton jets in each decay chain. There is a clear edge in 
their invariant mass distribution at M Xo , as shown in Fig. 25. This provides an absolute 
mass measurement and helps to remove some of the degeneracies discussed in the literature 
[50]. In addition, such reconstructions can be very useful in other precision measurements 
of the properties of the MSSM superpartners. For example, since we now have information 
about the LSP mass, and the direction of its decay products, it is easier to reconstruct the 
kinematics of the full event. In fact, we can fully recover the kinematics of the event using 
the same reconstruction method mentioned in the case of Z°+jets associated production in 
section 6.1.5. Gaining such information will significantly improve the prospect of measuring 
the spin of the LSP, which can be very challenging in the conventional scenario. 

7. Conclusions 

In this paper we have explored numerous aspects of model building and collider phenomenol- 
ogy for models of dark matter in which a non-abelian dark gauge symmetry is spontaneously 
broken at a GeV. Assuming a minimal dark gauge group, Gdark = SU(2) x U(l), we have 
surveyed a broad class of dark Higgs sectors that break enough symmetries (charge and cus- 
todial) to be phenomenologically viable. Furthermore, in order to accommodate the XDM 
and/or iDM scenarios, we have included the DM as the lightest state of a gauge multiplet of 
Crdark; which we have shown can naturally acquire mass splittings of ~ 100 keV — 1 MeV, 
from radiative corrections and higher-dimensional operators. 

We have also argued that an attractive option is to include supersymmetry, since this 
automatically generates a GeV symmetry breaking scale for the dark sector. In particular, 
we have proposed a novel mechanism whereby the U{1) kinetic mixing alone generates an 
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effective FI term for the dark hypercharge D-term, which in turn can induce SUSY breaking 
or a super-Higgs mechanism in the dark sector at a GeV. We emphasize that these GeV 
scale contributions will be always present, even if we choose to communicate SUSY breaking 
to the dark sector via more conventional gauge or gravity mediation. Also, in this scenario 
the dark matter multiplet can be uncharged under the SM gauge group, allowing us to evade 
astrophysical constraints which may arise if the dark matter can annihilate directly to 
(which in turn decay to p and 7r°). 

We have also given a detailed analysis of the gauge mediation scenario in which the 
dark matter is a 5 + 5 that communicates SUSY breaking from the SM to the dark sector, as 
originally proposed in Ref. [28]. Several benchmarks models with and without supersymmetry 
have been explicitly constructed here as examples. 

Our model building effort is far from exhaustive. Many directions remain to be explored. 
For example, in the context of supersymmetric implementations, we have assumed a natural 
solution for the n/Bfx problems in the dark sector. More detailed model building is certainly 
necessary to explore this issue further. 

We have also presented a broad analysis of possible collider signatures of a non-abelian 
dark sector. These results are completely independent of the specific details of our benchmark 
models (for example how SUSY is broken, or the exact choice of dark gauge group). Indeed, 
our assumptions involved only the existence of a kinetic mixing and the non-abelian gauge 
symmetry of the dark sector. 

We have found that dark gauge bosons can be produced via kinetic mixing with the 
photon in processes analogous to prompt photon production in the SM (with an additional 
suppression of e 2 on the rate). At the same time, this kinetic mixing also implies that dark 
gauge bosons can be produced by rare Z decays with a branching ratio of the order e 2 . The 
dark sector states dominantly decay into SM e ± and / u ± . With GeV invariant masses and 
typically large boosts of ~ 10s GeV, these leptons typically form highly collimated "lepton 
jets" . Indeed, a key feature of a non-abelian dark sector is that it produces lepton jets with 
multiple leptons, easily > 2 and possibly even 8. We expect that identifying more than 
2 leptons in a given lepton jet will significantly suppress backgrounds and enhance signal 
observability. With this in mind, there is the promising possibility that such signals may 
be observed at the LHC or possibly even at the Tevatron. Prompt photon-like production 
has a larger rate than rare Z-decay, but the latter is a complementary process that provides 
more handles on the dark sector. The 2-body decays of the dark sector into SM leptons 
do not generally lead to displaced vertices. However, such displaced vertices are in general 
possible in the case of a 3-body decay. If the 3-body decay is further suppressed by mass 
splitting or if the lightest scalar is the lightest dark state, then such a decay leads to $t which 
is collimated with the lepton jet. Several such displaced vertices will produce uncorrelated 
impact parameters of decay products. 

In supersymmetric implementations, the lightest superpartner (LSP) of the MSSM will 
decay further into the dark sector states which inevitably leads to lepton jets. This offers 
the possibility of detecting direct electroweak gaugino production with enhanced reach both 
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at the Tevatron and at the LHC. We can also take advantage of the unique kinematics of 
such SUSY cascades to perform absolute mass measurements of the MSSM's LSP. Missing 
energy is in general present in these events even if the gravitino is lighter than the lightest 
dark superpartner (LDSP), since the decay of the LDSP^ G + X has a lifetime much longer 
than the detector timescale. 

We would like to emphasize that we have only outlined the leading features of the collider 
phenomenology of this scenario. More detailed simulations are needed to obtain fully accu- 
rate estimates of the SM background with multiple collimated leptons. There are of course 
variations of the models considered here which could have new types of signals. An R-parity 
violating scenario could have more intriguing decay chains into the dark sector [28] . The small 
coupling between the observable sector and the dark sector could come with other types of 
"portals", such as the ones presented in Refs. [1, 51], which could lead to new collider signals. 
We observed that the modification of the SUSY decay chain will in general lead to better 
kinematical reconstruction. Such new information could result in precision measurements of 
other properties of the MSSM superparters involved in the decay chain. New techniques for 
taking advantage of such new information are certainly worthy of development. 
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A. Kinetic Mixing 

In this appendix we give a detailed description of kinetic mixing and its effect on dark sec- 
tor/SM couplings. To begin, we consider the non-SUSY case. Following the proposal of [1], 
we couple the dark sector to the SM via a gauge kinetic mixing between the dark and SM 
hypercharges (see Eq. (2.1)), much like what happens in the SM between the photon and the 
rho meson. This scheme is attractive because it does not break any symmetries of the SM 
and is hence less phenomenologically constrained. Moreover, since this operator is marginal, 
it can be generated at a very high scale, and will persist in the infrared. This implies that 
if both t/(l)'s are fundamental, the kinetic mixing is a UV boundary condition sensitive to 
physics at the highest scales. But if either U (1) is ultimately embedded in a GUT, kinetic 
mixing is only induced below the GUT scale by fields charged under both £7(l)'s. For exam- 
ple, by integrating out a multiplet of heavy fields <3?j of mass Mj that is charged under both 
dark and SM hypercharge, we find that 




(A-l) 
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where Qi and qi are the charges of <3?j under dark and SM hypercharge, and [i is the renor- 
malization scale. If for example, SM hypercharge is generated by symmetry breaking of some 
GUT group under which <E>j is charged, then J2i Qi = 0- ^ this multiplet has uniform qi 
charge, then the \x dependence cancels and the argument of the log becomes some ratio of 
scales in the multiplet, M/M' . For reasonable sizes of gy and g y , and a log contribution 
log M/M' ~ 1, this implies e ~ 1CT 4 - 1(T 3 , which is in the right range to explain DAMA. 

Next, let us consider how the kinetic mixing induces couplings between the dark sector 
and SM. At the electroweak scale, the terms involving the kinetic mixing are 

£gaugc mix = -\W 3IU/ Wf ~ \B^B^ - hb^f" + "-B^ (A-2) 

= -\z^ v - -F^ - hb^V" + ^(co S e w F^ - sin^V)^ (A-3) 

where F^ v and Z^ u are the fields strengths for the SM photon and Z boson, and in the second 
line we have gone from gauge eigenstate to mass eigenstate. Performing a field redefinition 
on the photon and the dark hypercharge gauge boson 

Afj! = Ay, - e cos 9 w by (A-4) 
b,/ = b lt + eamO w Zp (A-5) 

removes the kinetic mixing between the photon and Z, and removes the kinetic mixing be- 
tween the b and Z up to order e 3 . In addition, these shifts will modify the gauge-current 
couplings, AyJem + Z^Jg + byJ^ + w^Jw, as well as the gauge boson mass matrix. Since the 
photon is exactly massless, the shift of A has no effect on the mass matrix and simply couples 
b to the electromagnetic current of the SM. This is precisely the channel that will generate 
the leptons seen in astrophysical data. 

Analogously, the shift of b induces a coupling of the Z boson to the dark sector b current. 
However, unlike the photon, the b actually acquires a mass at ~ GeV, and furthermore mixes 
maximally with the dark iu's. For this reason shifting b induces a new mass mixing term 
between all the dark gauge bosons and Z^ of order em^/m| = (1 GeV/100 GeV) 2 x e. This 
in turn generates a mass suppressed coupling between b and the Z current and between Z 
and the w current! Thus, after removing the kinetic mixing, all the terms that couple the SM 
to the dark sector are 

-'-coupling = 

eby {cos9 w J? m + <D(ml/m? z )J%) + eZ^ (- sm6 w j£ + 0(m 2 b /m 2 z )J%) (A-6) 

where we have suppressed the mixing angles corresponding to the higher order contributions. 

If we now add SUSY, then the kinetic mixing becomes the expression shown in Eq. (4.1). 
This induces a mixing term for gauginos and D-terms. Since we have already considered the 
D-term mixing in Section 4.1, we focus on the gauginos. The new term is 

£ g augino mix = -ZieX^dy + h.c. (A-7) 
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where A^ and A^ are the dark and MSSM bino, respectively. Once again, since the dark bino 
is effectively massless at the electroweak scale, it is natural to shift it by A^ — > A^ + eA^. This 
yields a coupling term 

^coupling = e (\jjJb + 0(M~ b /Mz)\~ b JB) (A-8) 

i 

Jb = 9YY,QiH\Hi (A-10) 

i 

where J;, and Jb are the fermionic components of the dark and SM hypcrcharge supercurrents. 
Here the term that is O(M^fM^) arises from new mass mixing terms that arise from the 
gaugino shift. 



B. Conditions for Charge Breaking 

In general it is straightforward to achieve "electroweak" breaking for Gdark = SU{2) x U(l), 
simply by introducing a negative mass squared at the origin of Higgs field space. In the 
non-SUSY case this tachyon is inserted by hand, while in the SUSY case it arises naturally 
in the gauge mediation or kinetic mixing mediation scenarios mentioned in this paper. 

However, breaking G^ark completely, i.e. breaking charge, is a more difficult task. To see 
this, let us first consider the two Higgs doublet model. We can parameterize the vevs by 

h 1 =v 1 ( CMa ), h 2 = vJ°) (B-1) 
\ sin a I V / 

where h\ and hi are 2_ x / 2 and 2^/2-, respectively. For simplicity we have assumed that CP 
is preserved, and we have applied an SU(2) transformation to rotate hi into a single real 
component. For a given Higgs potential it is possible to determine whether charge is broken 
by considering the effective potential for the charge breaking angle a. Charge is preserved 
only if a = or tt at the minimum of the potential. Since \h\\ 2 and I/12I 2 are independent 
of a, only two renormalizable potential terms can contribute: \h\eh2\ 2 and hjehi- Naively, 
|/i|/i2| 2 contributes as well, but this term can be written as |/ij/i2| 2 = |^i| 2 |^-2| 2 — \h\eh^\' ■ 
Expanding these contributions in terms of a, the effective potential becomes 

V cS (a) = -^Acos 2 a + Bcosa (B-2) 

where A and B are a function of vip and the couplings. There is an extremum at a = 
arccosi?/A Checking that this point is stable, we find that a necessary condition for charge 
breaking is A < and \B/A\ < 1. 

Next, let us consider this in an example. In the MSSM, the quartic couplings are fixed 
by the D-terms, which turns out to fix A = g 2 > 0, which is why charge is left unbroken. 
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However, it is possible to push A below zero by introducing appropriate quartic contributions 
to the MSSM. In our SUSY benchmarks we accomplish this by including triplets. The one 
triplet SUSY benchmark has the superpotential: 

W = /x$Tr + fiH?eH 2 + \Hfe<S>H 2 (B-3) 

Let us consider the case where the triplet is heavy and we can integrate it out; this yields 
an effective theory of doublets in which charge breaking is simply determined. Our SUSY 
benchmarks are not in this decoupling limit, but nonetheless the physics of charge breaking 
in the low-energy doublet model appears to persist even as the triplet mass is lowered. In- 
tegrating out the triplet yields a quartic for the doublet with a coupling of A 2 m|/^|, where 
to$ is the soft mass for the triplet. Thus in V e g(a) for this model we find 

A = 5 2 + A 2 ^f (B-4) 

We see that A < only if there is a negative soft mass for the triplet that is appropriately 
large. While this can be easily engineered using gauge mediation, this not always possible 
generically. For example, if SUSY breaking is communicated in this theory via kinetic mixing 
mediation, then m| is not generated at the leading order, since <I> is a singlet under the 
dark hypercharge. On the other hand, we can easily remedy this by charging <I> under dark 
hypercharge — however, for anomaly cancellation we must also introduce a second triplet of 
opposite charge. In this model of two complex triplets, the kinetic mixing mediation will 
generate soft masses for the triplets. 

C. Z2 Symmetry in the tan/? = 1 Limit 

When tan/3 = 1, there is an enhanced 7L 2 symmetry of the two Higgs doublet model that 
makes it phenomenologically inviable. In particular, under this symmetry the w±„ have 
charge —1 and the and (dark photon) have charge +1. Since only z and a contain a 
component of the dark hypercharge, 6, only these states couple to SM electric charge. On the 
other hand, transitions between the different dark matter states are mediated by w± alone. 
Thus it is necessary to break this Z 2 . 

The origin of this % 2 is as follows. When tan/? = 1, then h\ and h 2 have the same 
magnitude; thus, they can be simply thought of as two spinors that correspond to two unit 
vectors in 3-space. Next, h\ and h 2 uniquely define a third direction which bisects the angle 
a between them. Rotations of 180° around this axis, followed by h± h 2 , leave the vacuum 
invariant. Thus, all states in the low-energy theory are eigenstates of this Z 2 . Since this Z 2 is 
a subgroup of SU(2), it acts nicely on {wi, w 2 , W3}. It is obvious by choosing a basis where 
W3 points along the axis of rotation, that two of the SU (2) gauge bosons are odd under this 
7j 2 and the remaining one is even. Thus, the latter is the only state that can mix with the 
dark hypercharge, b, since it is also Z 2 neutral. 
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D. Supersymmetric Contribution to Mass Splitting 



In this appendix we present the supersymmetric contributions to the mass corrections of a 
Dirac fermion, We show that in the parameter region we are interested in, those are 
negligible. These results are well-known (see for example [34]) and are presented here for 
completeness. We begin with the non-supersymmetric contributions. Consider, therefore, a 
theory with a SM-like weak gauge-group SU(2) x U(l) which in general is broken down to 
nothing at some scale. Also, for simplicity, we take the Dirac fermion to be charged as 2]/ 2 
with a mass My much larger than the Higgs scale of the theory. Similar conclusions hold for 
any representation of the gauge-group. At low energies, the masses of the two components 
are split. There are both wave-functions and mass insertion diagrams, given by, 

A 

= - % ^f f dxxlogA (D-l) 





X " 



^M 9 Jdx\ogA (D-2) 



where A = ((1 — x) 2 M| + xM^j , M A is the mass of the gauge boson, and gA is its coupling 
to the fermion. We neglected the divergent part since it cancels when considering the mass 
splitting. The propagating gauge boson is any one of the four massive vector bosons. 

For the rest of this section we consider the simple case where one gauge boson is left 
massless. In that case we have, 

aM 9 f f xMl \ m*>m z aM z 
5M 9 = —— dx\og[l + - -ttth) >^r~ ( D - 3 ) 



2tt J b V (1 " x) 2 Ml 

Where a is associated with the massless gauge boson coupling to matter. For a general 
multiplet of the gauge group the mass splitting between two eigenstates, i and j of T 3 is given 
by, 

5M tJ = ^( q f- q ])M z (D-4) 
- y (ffi 3 ) 2 " (Tf) 2 ) (M z - M w ) , 

where the notation is explained after Eq. (3.3) in Section 3. When the "photon" is also 
massive the correction goes as the splitting between the gauge boson masses, 

oi 

5My « - (M z - M 7 ) + . . . (D-5) 

However, the precise formula requires the vector boson mass eigenstates and is not simple in 
general. 
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The supersymmetric contribution is through a similar loop to the wave-function renor- 
malization above only with a gaugino - slepton loop replacing the gauge boson - lepton 



propagators 15 , 




ie 2 
8^ 



dx(l - x) log A (D-6) 



with A = ^(1 — x) 2 M? + xMj^J . The contribution to the splitting is then, 

Q M, f ( (1 - x) 2 M 2 -, + xM 2 ~ \ 

SM^ = ^ dxlog V - Lr-%± -J (D-7) 

2vr J x) 2 M| + xM? J 

where we used M™, My ( M^o) casually to designate the charged and neutral gauge 

bosons (leptons). Clearly, in the limit where all the masses are equal the integral vanishes. 
Therefore, the only contribution to the splitting comes from possible differences between M^ 
and M~ (or and M^ ) . If we denote the soft supersymmetry contribution to the gaugino 
masses by M\ we have, 

M^-M?^ 2M X M W (D-8) 

Therefore, this contribution to the mass splitting in Eq. (D-7) is suppressed by M\/M 9 
compared with the non-supersymmetric contribution in Eq. (D-5). 

Another possibility is that the charged slepton is split from the neutral one by SU (2) 
D-term contributions. However, those contributions are to the mass squared. Writing — 
Mq w Myy we see that again the contribution to the mass splitting is suppressed by Mw/M 9 
with respect to Eq. (D-5). 
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